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Abstract. We discuss the role of quantum chaos in atomic nuclei. After 
reviewing the basic assumptions of the nuclear shell model, we analyze the 
spectral statistics of the energy levels obtained with realistic shell-model 
calculations in the fp shell. In particular, for Ca isotopes we observe a 
transition from order to chaos by increasing the excitation energy and a 
clear quantum signature of the breaking of integrability by changing the 
single-particle spacings. 



1 Chaos and Quantum Chaos 

Quantum chaos is the study of the properties of quantal systems which are 
classically chaotic, thus with exponential divergence of initially closed tra- 
jectories in the phase space [1]. 

The energy fluctuation properties of systems with underlying classical 
chaotic behaviour and time-reversal symmetry agree with the predictions of 
the Gaussian Orthogonal Ensemble (GOE) of random matrix theory, whereas 
quantum analogs of classically integrable systems display the characteristics 
of the Poisson statistics [2]. 
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The most used quantity to study the local fluctuations of the energy 
levels is the spectral statistics P(s). P(s) is the distribution of nearest- 
neighbour spacings Sj = (E i+1 — Ej) of the unfolded levels Ei. It is obtained 
by accumulating the number of spacings that lie within the bin (s, s+ As) and 
then normalizing P(s) to unity. For quantum systems whose classical analogs 
are integrable, P(s) is expected to follow the Poisson limit, i.e. P(s) = 
exp (— s). On the other hand, quantal analogs of chaotic systems exhibit 
the spectral properties of GOE with P(s) = (7r/2)sexp (— f s 2 ), which is the 
so-called Wigner distribution. The distribution P(s) is the best spectral 
statistics to analyze shorter series of energy levels and intermediate regions 
between order and chaos [1,2]. 



One of the best systems for the study of quantum chaos is the atomic nucleus. 
In fact, its experimental energy levels have been studied in the domain of 
neutron and proton resonances near the nucleon emission threshold, where a 
large number of levels with the same spin and parity in the same nucleus are 
present, and an excellent agreement with GOE predictions has been found 



In the nuclear shell model [4], the nuclear states, like the electronic states 
in atoms, are described in terms of the motion of nucleons in a mean-field. 
But, while the nuclear field is generated by the interactions of the nucleons 
(protons and neutrons), the atomic field is mainly governed by the interaction 
of the electrons with the nucleus. 

The Hamiltonian of the shell model can be written as: 



where A is the number of nucleons, rrii is the mass of the i-th nucleon, and 
Ui is the mean-field. The choice of the mean-field is crucial; C/j may be 
obtained by the usual methods using the many body theory [4]. 

To obtain a good agreement between the shell model results and the 
experimental data, it is necessary to add a residual interaction Hr so that 
the total hamiltonian H can be written: 



2 The Nuclear Shell Model 



[3]. 




(1) 



H = if o + Hr , 



(2) 
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where Hr is the part of nucleon-nucleon interaction not included in Hq. 
Using second quantization formalism [4], we can write: 



#o = E 




(3) 



Q 



(4) 



where the labels denote the accessible single-particle states, e a is the corre- 
sponding single-particle energy, and < a/3|y|57 > is the two-body matrix 
element of the nuclear residual interaction. The operators a+ and a a are 
the fermionic creation and annihilation operators of the ct-th single nucleon 
state and Hr can be calculated by the Hartree-Fock equations, starting from 
the free nucleon interactions. Many nucleons are frozen in the deeper shells 
of the mean field potential and form an inert core; only a few nucleons par- 
tially populate the single particle shells outside the core. These are called 
valence-nucleons. So there are N valence-nucleons, m active shells and a 
finite number of energy levels. It is standard procedure to cut the basis 
states by introducing a finite number of configurations which are sufficient 
to describe the first excitation states [4]. 



3 Spectral Statistics and the Brody Distribu- 



We study the (f7/2,P3/2,f5/2,Pi/2) shell-model space, assuming 40 Ca as an 
inert core with 4 < N < 10 valence-nucleons. The diagonalizations are 
performed by using Lanczos and Householder algorithms with the code AN- 
TOINE [5,6]. For a fixed number of valence protons and neutrons we cal- 
culate the energy spectrum for total angular momentum J and total isospin 
T. The single-particle energies (in MeV) are e 7 / 2 = 0, e 3 / 2 = 2, ei/ 2 = 4 
and e 5 / 2 = 6.5. The residual interaction we use is a minimally modified 
Kuo-Brown realistic force with monopole improvements [6]. 

For the low-lying levels, the spectrum is mapped onto unfolded levels 
with quasi-uniform level density by using both the constant temperature 
formula [7] and the local unfolding method [8]. The two procedures give 
the same results. For the full spectrum we use the local unfolding because 



tion 
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the constant temperature formula is valid only when the level density grows 
exponentially. 

To quantify the chaoticity of the distribution P(s) (nearest-neighbour 
spacings of the energy levels) in terms of a parameter, we compare it to the 
Brody distribution 

P(s,u) = a(cu + l)s w exp(-as w+1 ) , (5) 

with 

a = (r^ir 1 . (6) 

This distribution interpolates between the Poisson distribution {uj = 0) of 
regular systems and the Wigner distribution (u — 1) of chaotic ones (GOE). 
The parameter uj can be used as a simple quantitative measure of the degree 
of chaoticity [9]. 

Table 1. Brody parameter u for Ca isotopes. 



44 Ca 


45 Ca 


46 Ca 


47 Ca 


48 Ca 


49 Ca 50 Ca 


0.69 


0.75 


0.99 


0.98 


0.95 


1.00 0.87 



Table 1 shows the Brody parameter uj for the whole spectrum of the 
analyzed Ca isoptopes (only neutrons in the fp shell), which range from 
44 Ca to 50 Ca. We see that only the lightest Ca isotopes are not fully chaotic. 

It now becomes interesting to analyze the Brody parameter as a function 
of the excitation energy. We calculate the P{s) distribution and the Brody 
parameter up to a fixed value of the excitation energy above the yrast lines. 
Obviously this cannot be done for the lightest Ca isotopes because few levels 
are involved. 

Table 2. Brody parameter wasa function the the Energy for Ca isotopes. 



Energy (MeV) 


48 Ca 


49 Ca 


50 Ca 


6 


0.54 


0.63 


0.72 


8 


0.65 


0.70 


0.79 


10 


0.78 


0.81 


0.84 


12 


0.83 


0.82 


0.86 


14 


0.90 


0.86 


0.88 


16 


0.87 


0.92 


0.90 


18 


0.93 


0.95 


0.92 
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The results are written in Table 2 for 48 Ca, 49 Ca and 50 Ca and show a 
strong energy dependence: there is an order-chaos transition by increasing 
the excitation energy. 

Another important aspect is the effect of the one-body Hamiltonian 
on the two-body residual interaction when the single-particle spacings are 
changed. 



Table 3. Brody parameter u for Ca isotopes. 



Single Particle Spacings 


44 Ca 


45 Ca 


degenerate 


0.85 


0.98 


normal 


0.69 


0.75 


double-spaced 


0.29 


0.58 



Table 3 shows that with degenerate single-particle levels the isotopes Ca 
and 45 Ca are chaotic while for double-spaced single-particle levels they are 
quasi-regular. This interesting effect is a clear quantum signature of the 
breaking of the integrability due to the residual interaction. In fact, the one- 
body Hamiltonian is classically integrable because it is the sum of harmonic 
oscillators while the two-body residual interaction is strongly non-linear. By 
increasing the single-particle spacings, the single-particle mean-field motion 
in the valence orbits suffers less disturbance and is thus more regular. 

4 Summary 

We have studied quantum chaos in atomic nuclei by using the nuclear shell 
model. We have calculated the spectral statistics of Ca isotopes by using 
Lanczos and Householder algorithms with the shell-model code ANTOINE. 
The Brody parameter, which fits the nearest-neighbour spacings distribution 
of the energy levels, shows an order-chaos transition when the excitation 
energy is increased and when the single-particle spacings are changed. 
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